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Singly charged, positive and negative, atomic ions have been observed to sublime freely from 4d 
and Sd transition metal surfaces in the 1800 to 2600 K temperature range. The rates of sublimation 
were found to be about equal to those which would be observed under conditions of thermal equilib- 
rium. Ion sublimation energies were obtained from the temperature dependence of these sublimation 
rates. The energies were found to be consistent with currently accepted values for the atom sub- 
limation energy, electron work function, and first ionization potentials of the elements investigated. 

The electron affinities of these metal, atoms were determined from an isothermal measurement 
of the ratio of positive to negative ion sublimation rates. The results were compared with some recent 
theoretical estimates of this quantity for those Sd transition metals whose valence electron configura- 
tions are analogous to the ^d and 5r/ elements considered here. 
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An atomically clean metal surface may be char- 
acterized by the energies associated with the sublima- 
tion processes: 

(1) e 8 — > e g ; (/>, the electron work function 

(2) M ii —>M !f ; /o, the atom sublimation energy 

(3) M, s — e s -^ M*j,; / + , the positive ion sublimation 

energy 

(4) M s -\-e s —> Mj ; ; /_, the negative ion sublimation 

energy. 

At temperatures less than about § of the melting 
point, the energies associated with these processes 
may be strongly influenced by the geometric arrange- 
ment of the surface atoms. The observed difference 
in work function between different crystalline faces 
is an example of such anisotropic behavior. At high 
temperatures the geometric structure and symmetry 
characteristics of the surface become less important 
as a result of rapid surface migration and large sub- 
limation rates. In fact, in the neighborhood of the 
melting point the surface can be thought of as a two- 
dimensional liquid whose properties tend to be uni- 
form in contrast to the considerable anisotropy 
exhibited at lower temperatures. 

Process (3) has been observed for a number of 
refractory transition metals [1-8, 10]. i Early attempts 
at a reliable measurement of /+ were limited by primi- 



1 Figures in brackets indicate the literature references at the end of this paper. 



tive vacuum conditions and the presence of easily 
ionized alkali atom impurities. Process (4) can be 
observed when A(M), the electron affinity of A/, is 
positive {Mg is stable). It has been reported for 
graphite, tungsten, rhenium, molybdenum, and tan- 
talum [7, 8, 9, 10]. 



1. A Generalized Saha-Langmuir Equation 

Consider the thermodynamic equilibrium existing 
inside of an isothermal metal enclosure. In addition 
to blackbody photons, the stable gaseous species 
are electrons e fh metal atoms M lh singly charged 
positive ions Mt } •> and singly charged negative ions My 
when the electron affinity of the atom is positive. 
If the number density of charged particles is suffi- 
ciently small so that space-charge effects may be 
neglected, then the partial pressures (p x ) are governed 
by the following equilibria: 

My + e (J ^ Mf; AE = + I upper sign for positive 

ionization 
— — A lower sign for negative 
ionization. 

The quantity A£ is the energy (or free-energy) change 
for these ionization processes at absolute zero. The 
equilibrium constants are given by 

(p Jpo) (Pe) ±x = (QJQo) (Qe)^exp (- A£/AT) . (1) 
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The Q's are the molar partition functions for these 
gaseous species which are assumed to be in thermal 
equilibrium with the walls of the enclosure. The 
partial pressure of electrons (p e ) is given by the 
well-known Richardson equation: 



Pe = Qe exp (—<j>lkT) 



(2) 



The electronic work function <j> is defined as the 
free-energy change at K for electron sublimation. 
Substituting eq (2) into (1) yields the result that 

(p±lpo) = (Q±IQo) exp [{-t±E±4>)lkT]. (3) 

The flux (j) of particles evaporating from (or con- 
densing upon) the walls of the enclosure is given by 



j=apl(27rmkT) 1 / 2 



(4) 



where a is the evaporation (or condensation) coeffi- 
cient and p is the equilibrium vapor pressure. Since 
the masses (m) of the atomic ions and atoms are 
very nearly equal and, assuming that the evaporation 
coefficients for ions and atoms are also equal, the 
ratio of sublimation rates obtained from eqs (3) and 
(4) is 



j±lj =(OJQo) exp [(-\E±4>)lkT]. 



(5) 



For the atomic and ionic species under consideration, 
the ratio of molar partition functions is equal to the 
ratio of electronic partition functions: 

«?±/<?o) = («±/«o), 

since the translational factors cancel because of the 
equality of the particle masses. When Russel-Saunders 
spin-orbit coupling is obeyed, the electronic partition 
function can be written as 



(o = ^g n exp (—eJkT) 



(6) 



where e = and g n = 2J n + l. The quantity g n is the 
statistical weight of the level of energy e n , and J n the 
inner quantum number of the electronic state n. 
The summation need only be carried out for those 
states which are significantly populated at the tem- 
perature T. Under the condition of thermal equilibrium 
on the metal surface, eq (5) yields a result for the subli- 
mation ratio of singly charged positive ions to neutral 
atoms which may be recognized as the familiar 
Saha-Langmuir equation [11]; namely, 



j+ljo=(a)+l<Oo) exp [((/> — I)/kT]. 



(7) 



The ratio of positive to negative self-surface ionization, 
obtained from eqs (7) and (8), is 

j+//- = (»+/©-) exp [(24>-I-A)lkr\. (9) 

The atomic properties / and A can be shown to be 
related to the sublimation energies (/>, / , /+, and /_ 
from a consideration of the cycles given in figure 1; 
namely, that 



0-/=/o-/ + , 



A — (j> = l — /_, 



and 



20-7-/4 = /_-/+ = / + (/>-/*-/+. 



(10) 



(11) 



(12) 



The equilibrium vapor pressure of metal atoms (p ) 
is given by 



Po = exp [(AS/ft)-(/ /AT)L 



(13) 



where AS and / are the atom sublimation entropy and 
energy at some temperature T. Substituting eq (13) 
into (3), one obtains for the sublimation pressure of 
positive and negative ions: 



P± 



Qo 
Q± 



exp 



-AS 



exp 



-Ag±(/>-/o 
kT 



(14) 



Using eq (4) and collecting terms, the dependence of 
the rates of positive- and negative-ion sublimation 
upon temperature is given by 



""©-(-¥) 



— exp l-j^- 



(2tt mk) 



(15a) 




METAL SURFACE 





Similarly, when negative self-surface ionization is 
considered, eq (5) yields the ratio of subliming negative 
ions to atoms: 



j-ljo= (o)-/o) ) exp [(A — <l>)lkT]. 



(8) 



PSI 



NSI: 



JL« 



^o 



FIGURE 1. Thermal cycles for positive and negative self-surface 
ionization. 
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and 



J- 



7*1/2 



exp , — 



(277mA-) 1 ' 2 



exp 



A-<t>-l u 
kT 



(15b) 



where from eqs (10) and (11) it is seen that 
<f) — I — [{) = — 1+ and A — <£> — Zo = — /_. The /'s are sub- 
limation energies at some mean temperature T in 
the middle of an experimental temperature range, 
and, like (/>, /, and A, remain virtually unchanged from 
their values at absolute zero. 

2. Apparatus 

The apparatus used in these studies consisted of a 
direction-focusing mass spectrometer. The surface 
ionization source of the mass analyzer was designed 
to facilitate the collection of negative ions, which 
evaporate in observable quantities only at surface 
temperatures at which electron emission is intense. 
A cylindrical modification of a Pierce lens [12] was 
used to collect the charged particles evaporating from 
the experimental surfaces. When the apparatus was 
biased for the detection of negative charges, the large 
electron flux focused by the lens was deflected by a 
small magnetic field and collected on a molybdenum 
plate. Because of space limitations in the ion source, 
the lens design was limited to the efficient collection 
of only about 2 mA of electron emission. At surface 
temperatures in excess of that which yields this value 
of the saturation electron current, a negative space 
charge begins to accumulate and the collection effi- 
ciency for negative ions decreases. This problem is 
not encountered in the case of positive ion sublimation 
so that meaningful measurements at a constant ion 
collection efficiency can be made at higher surface 
temperatures and larger positive ion currents. 

Both the positive and negative ions were accelerated 
by the components of the Pierce lens to an energy of 
4 kV prior to their traversal through the entrance slit 
of the 15-cm radius, 60° sector, magnetic analyzer. 
The first dynode of a 16-stage Cu-Be electron multiplier 
detector was maintained at ground potential, so that 
both positive and negative mass analyzed ions were col- 
lected with 4 kV of energy. The multiplier was operated 
with 190 V between stages so that the sixteenth dynode 
was at +3040 V for both positive- and negative-ion 
detection. The output pulses were 20 to 30 ns. wide and 
were amplified with a voltage gain of 5000. The re- 
sultant pulse-height distribution was found to be about 
the same for both positive and negative ions. Thus, the 
secondary electron yields for 4 kV positive and nega- 
tive ions impacting on the first dynode surface do not 
appear to exhibit any significant charge discrimination. 
This is of particular importance for the determination 
of the ratio of positive to negative ion sublimation rates. 
The amplified pulses were then processed to give 
uniform 1-V, 50-ns pulses which actuated a counter 
with a 10-MHz bandpass. The background count for 



this instrument was usually less than 10 pulses/minute 
with the magnetic analyzer detuned and the experi- 
mental surface temperature well below 1800 K. 



3. Experimental Results 

The temperature dependence of the rates of sub- 
limation of singly charged positive and negative ions 
of tungsten, rhenium, molybdenum, and tantalum are 
shown in figures 2, 3, and 4. Figure 5 gives the results 
of some measurements of the rate of positive ion sub- 
limation from a niobium surface. The quantity j± 
(corrected) is the experimental ion count corrected for 
the slightly temperature dependent factor 

[7 T1 / 2 (o) /co ± ) exp (-AS/A)] 

given on the lefthand side of eqs (15a) and (15b). The 
values obtained for this factor were normalized in the 
middle of each experimental temperature range. The 
values for AS were obtained from the compilation of 
Stull and Sinke [13J. The o>o and co+ were calculated 
from eq (6) using Moore's energy level compilation 
[14]. For the negative ion electronic partition function 
o>_, only one bound state was assumed with a./o value 
equal to that of the ground state of the atom with the 
next highest atomic number. This assumption appears 
to be reasonable since virtually all of the negative ions 
which have been observed so far exhibit this property 
[15]. The evaporation coefficients (a±) were assumed 
to be temperature independent with values of unity. 
Table 1 gives a summary of the sublimation energies, 
obtained in this laboratory, for the positive ions of a 
number of transition metals [5, 7, 8, 10] of rows 4 and 
5 of the periodic table. The agreement between the 
/+ calculated from eq (10) and the experimental values 
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FlGURE 2. Positive and negative self-surface ionization of tungsten 
and rhenium. 

The uncertainties in the sublimation energies are the standard errors in the slopes of 
the lines. 
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FIGURE 3. Positive and negative self-surface ionization of 
molybdenum. 

The error bars represent the range of values obtained in at least ten determinations. The 
uncertainty in the sublimation energy is the standard error in the slope. 



is seen to be within the estimated uncertainties. The 
quantities /o, /, and (/) for these metals were obtained 
from recent compilations [10, 13, 14, 16]. The intercepts 
of all of the experimental lines given in figures 2 to 5, 
divided by the subliming surface areas, fall within a 
factor of three of a value 5 X 10 28 cm -2 sec 1 . This is 
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Figure 4. Positive and negative ion sublimation from a tantalum 
surface. 

The error bars represent the range of values obtained in at least ten determinations. The 
uncertainty in the sublimation energy is the standard error in the slope. 

consistent with a maximum limiting sublimation rate 
of one monolayer (~ 10 15 cm -2 ) per lattice vibration 
time (~ 10" 14 s). It can therefore be concluded that a 
relatively large fraction of the evaporating ions was 
focused by the Pierce lens into the mass analyzer and 
ultimately detected by the electron multiplier. This 
result, together with the demonstrated validity of the 
thermal cycle eq (10), supports the conclusion that 
the observed free sublimation rates are approximately 
equal to those which would be observed in a closed 

Table 1. Positive ion sublimation energies 



Ion 


T 


t+(obs.) 


b (/o + /-c/>) 


93 Nb + 

98 Mo + 
181 Ta + 

186 W + 

187 Re + 


(K) 
2100 
1925 
2280 
2450 
2130 


(eV) 
10.4 
9.5 
11.2 
12.1 
10.7 


(eV) 
7.5 + 6.8-3.9=10.5 

6.7 + 7.1-4.3 = 9.5 
7.9 + 7.9-4.1 = 11.7 
8.9 + 8.0-4.5=12.4 

7.8 + 7.9-5.0=10.7 



a Typical standard errors in the least square determination of /+ 
are between 0.1 and 0.2 eV. 

b Uncertainties in / , / and <f) are each about 0.1 eV. 
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THE EMISSION OF Nb FROM AN 
INCANDESCENT NIOBIUM SURFACE 



* 2I0 ° X 

v£ + = -( 10.39 ± 0.15) eV 




4.6 4.7 4.8 4.9 5.0 

f IK)"' 

FIGURE 5. The sublimation o/ 93 Nb + /rom a heated niobium 
surface. 

The uncertainty in the sublimation energy is the standard error in the slope of the line. 



Table 2. Isothermal positive to negative ion ratios obtained for 
the self-surface ionization of molybdenum, tantalium, tungsten, 
and rhenium. 

The errors given lory, and / are average deviations from an arithmetic mean of at least 

si\ successive ion counts. 



Ion 


T 


i+ 


;- 


AT In [co_y+/a> +./'-] 




(K) 


(ions/sec) 


{ions/ sec) 


(eV) 


98 Mo 


1903 


7.4 ±0.4 


0.30±0.12 


0.53 ±0.07 




1924 


13.3 ±0.6 


0.45 ±0.08 


0.56 ±0.03 




1945 


22.3 ±0.5 


1.40±0.20 


0.46 ±0.02 




1965 


37.5 ±0.4 


1.87±0.13 


0.51 ±0.01 




1985 


66.8 ±1.5 


2.30±0.15 


0.58 ±0.01 




0.53 ±0.03 eV 


181 Ta 


2000 


0.18 ±0.04 


0.03 ±0.01 


-0.04 ±0.08 




2047 


0.58±0.10 


0.10±0.03 


-0.06 ±0.06 




2077 


1.67±0.13 


0.18 ±0.05 


+ 0.02 ±0.06 




-0.03 ±0.07 eV 


18«W 


2225 


3.0±0.2 


0.07 ±0.03 


0.76 ±0.08 




2250 


6.6 ±0.2 


0.18±0.03 


0.74 ±0.04 




2280 


12.6±0.2 


0.30 ±0.06 


0.77 ±0.04 




0.76 ±0.05 eV 


.H7 Re 


2272 


7,210±40 


0.32 ±0.05 


2.01 ±0.03 




2290 


11,020 ±60 


0.49 ±0.04 


2.02 ±0.02 




2314 


18,650 ±100 


0.84 ±0.08 


2.04 ±0.02 




2333 


30,730 ±120 


1.40±0.15 


2.06 ±0.02 




2.03 ±0.02 eV 



system in thermal equilibrium with vaporization 
coefficients (a ± ) not very different from unity. 

The ratio of positive to negative self-surface ioniza- 
tion could be reliably measured only over a rather 
limited temperature range. This is seen in figures 2 to 
4, where the lowest temperature was determined by 
the smallest negative ion sublimation rate detectable 
with the present apparatus (5 to 20 ions/minute). The 
upper temperature was limited by an electron emission 
which was sufficient to saturate the Pierce lens. This 
resulted in a negative space charge buildup with a 
consequent lowering of the negative ion collection 
efficiency. It is seen that this limiting temperature was 
about 2100 K for the low work function metals, tan- 
talum and molybdenum, while the higher work func- 
tions of tungsten and rhenium allowed efficient nega- 
tive ion collection up to temperatures of 2300 and 2500 
K, respectively. 

The results obtained for these positive to negative 
ion ratios are given in table 2. The quantities j+ and 
j_ are the actual ion count rates observed at each of 
the experimental temperatures. The partition function 
ratios (co-/co + ) were calculated in the manner de- 
scribed above. The quantity k T In [(a>-/a> + ) (/+//-)]» 
given in the last column, is equal to (2<f>—I—A) or 
(0 + Zo — /+— ^) according to eqs (9) and (12). Solving 
(9), (10), (11), and (12) for A, the electron affinity, yields 



Table 3 gives the electron affinities for the four 
transition metals. They were calculated using eq (16). 
A reliable value for /_ could be obtained from eq (15b) 
only in the case of rhenium, as shown in figure 2. Con- 
sequently, it was possible to make two independent 
estimates of its electron affinity. It is seen that, de- 
spite the small absolute value of A (Re), the agreement 
between the two estimates is good and considerable 
confidence in the stability of Re" is gained. 

Table 3. The electron affinities of Mo, Ta, W, and Re as calculated 
from eq (16). 



A = <f> + lo—l- = <l> + lo-l+- 



-Win 



o>-J+ 

0)+j- 



(16) 



A {Mo): 
/f(Ta): 

A(W): 
A(Re): 



<f> + lo — l+ — kT\n [w-j+lcj+j-] 



4.3 + 6.7-9.5-0.5=1.0 
4.1 + 7.9-11.2 + 0.03 = 0.8 
4.5 + 8.9-12.1-0.8 = 0.5 
5.0 + 7.8-10.7-2.0 = 0.1 



(* + /„-/_) 



(eV) 



5.0+7.8-12.6 = 0.2 



4. Conclusions 

It has been shown that at temperatures in excess 
of about 1800 K the sublimation of singly charged, 
positive and negative ions from refractory metal sur- 
faces are in accordance with the thermal cycle eqs 
(10) and (11). The measured quantities /+, /_, and A 
are consistent with the currently accepted values for 



41 



/o, (/), and /. Under the high-temperature conditions 
of these measurements, no evidence for surface 
anisotropy was observed. 

A detailed review of atomic negative ions has been 
given recently by Smirnov [17]. He summarized both 
the experimental and theoretical electron affinity 
estimates that have been made for the elements of 
rows, 2, 3, and 4 and columns 1 and 6 of the periodic 
table. The experimental results given here are the 
first determinations of the electron affinities of the 
transition elements in rows 5 and 6. The method 
described should be applicable to all of those refrac- 
tory metals whose atoms can form stable negative ions. 

Clementi [18] has made a set of Hartree-Fock cal- 
culations for the 4th-row transition elements. From 
these calculations he was able to make correlation 
energy corrections to the 4th-row negative ion energies 
and hence obtained estimates of the electron affinities. 
Table 4 shows a comparison between these theoretical 
estimates and the self-surface ionization values de- 
termined in this laboratory for the 5th- and 6th-row 
transition metals. The value for niobium was estimated 
on the basis of approximate negative ion yields. Be- 
cause of its low work function [5] (~ 3.9 eV), copious 
electron emission made it impossible to obtain accu- 
rate positive to negative ion ratios with the present 
instrument. Molybdenum and chromium as well as 
tantalum and vanadium have analogous valence elec- 



Table 4. Electron affinities 
(Comparison with theoretical estimates) 



E. Clementi: 
Phys. Rev. 


4th row 
valence 
electron 
configur- 
ation 


V 


Cr 


Mn 


135, 

A980 
(1964). 
(correla- 
tion 
energy 
estimates) 


3dHs 2 


MHs 


3d 5 4s 2 




A(eV) 


0.94 ±0.25 


0.98 ±0.35 


-1.07 ±0.20 


This work: 
(self- 


5th row 
valence 
electron 
configur- 
tion 


Nb 


Mo 


Tc 


surface 
ioniza- 
tion) 


4d 4 5s 


4d 5 5s 


4d 5 5s 2 




A(eV) 


-1.0 (est.) 


1.0±0.2 










6th row 
valence 
electron 
configur- 
ation 


Ta 


W 


Re 




5d 3 6s 2 


5d 4 6s 2 


5d 5 6s 2 




A(eV) 


0.8±0.3 


0.5±0.3 


0.15±0.10 



tron configurations and there is remarkably good 
agreement between the theoretical and measured 
electron affinities. In the case of the manganese- 
rhenium pair, however, where there is a full s and 
half-filled d shell in their ground-state configurations, 
there is disagreement. Clementi predicts that such a 
valence electron structure should not yield a stable 
negative ion. In fact, his calculation gives a value of 
— 1.1 eV for ^(Mn). The negative self-surface ioniza- 
tion of rhenium has shown however, that the Sd 5 6s 2 
structure is indeed capable of producing a stable 
rhenium negative ion with a binding energy of about 
0.15 eV. 

A number of years ago, Gordy and Thomas [19] 
found that a linear relationship exists between the 
electronegativity and work function of metals. This 
empirical fact can be understood in terms of eqs (10) 
and (11), since elimination of Z yields 



= 



J+A L 



U 



(17) 



The quantity (I-\-A)/2 is the absolute electronega- 
tivity QC m ) first defined by Mulliken [20]. Gordy and 
Thomas showed that for most of the metals in the 
periodic table the following relation was obeyed: 



(/> = 2.3X p + 0.34, 



(18) 



where the quantity X p is the electronegativity given 
by Pauling [21]. Skinner and Pritchard [22] have 
shown that X m = 3A X p so that (18) may be written in 
terms of X m as follows: 



<$> = 0.74 X m + 0.34. 



(19) 



Comparing this result with eq (17) shows that, while 
the linear dependence of <j> on X m is predicted, the 
empirically determined coefficient is not quite unity 
as required by eq (17). However, detailed examination 
of the Gordy-Thomas correlation shows that a 25 per- 
cent increase in the coefficient of X m would not do 
much of an injustice to the data which they used. 
In any case, the results of the present work strongly 
suggests the general applicability of eq (17) rather than 
its empirical counterpart, eq (19). It should be pointed 
out that the "constant" term in eq (19) is really given 
by (Z_ — Z+)/2, which from eq (16) is seen to be equal 
to kT/2 In [aj-lj+/(o + j-]. As shown in table 2 this 
quantity depends upon the metal to which it is applied. 
Consequently, the empirical validity of eq (19) can only 
be considered approximate at best. 
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